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Results from the data obtained in the first physics run of the Relativistic Heavy Ion Collider 
(RHIC) have shown suprisingly large elliptic flow and suprisingly small HBT radii. Attempts to 
explain both results in a consistant picture have so far been unsuccessful. The present work shows 
that a simple thermal-like initial state model coupled to a hadronic rescattering calculation can 
explain reasonably well both elliptic flow and HBT results from RHIC. The calculation suggests a 
very early hadronization time of about 1 fm/c after the initial collision of the nuclei. 
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I. INTRODUCTION 

Results of the Year-1 running of the Relativistic Heavy 
Ion Collider (RHIC) for Au+Au collisions at y/s = 
130 GeV have shown suprisingly large pion elliptic flow 
101 and suprisingly small radii from two-pion Hanbury- 
Brown-Twiss interferometry (HBT) ||. Attempts to 
explain both results in a consistant picture have so far 
been unsuccessful. Hydrodynamical models agree with 
the large elliptic flow seen in the RHIC data M but sig- 
nificantly disagree with the experimental HBT radii f^. 
On the other hand, relativistic quantum molecular dy- 
namics calculations which include hadronic rescattering, 
for example RQMD v2.4 [Q, significantly underpredict 
the elliptic flow seen in the RHIC data |7] but predict 
pion HBT radii comparable to the data |^]. A calcula- 
tion has recently been made to extract HBT radii with 
a hydrodynamical model coupled with a hadronic rescat- 
tering afterburner with the result that the HBT radii are 
significanly larger than measurements j^. This lack of 
a single model to explain both results has been our first 
big mystery from RHIC. It has been suggested that we 
should call into question our current understanding of 
what information pion HBT measurements give us [[lO[ . 

In an effort to address this mystery, the present work 
explores a somewhat different picture of the nuclear 
collision than those presented above. In this picture, 
hadronization into a simple thermal-like state occurs 
soon (about 1 fm/c) after the initial collision of the nuclei 
followed by hadronic rescattering until freezeout. The 
goal will thus be to test whether hadronic rescattering 
alone can generate enough general flow in the system 
to explain both the elliptic flow and HBT results from 
RHIC starting from a simple non-flowing initial state. 
Note that this approach, while similar to the relativis- 
tic quantum molecular dynamics calculations mentioned 
above, differs from them in the choice of the model of 
the initial state of the system before hadronic rescatter- 



ing commences (e.g. the other models use a color string 
picture for the initial state) ||^. The simple initial state 
model used here is meant to be a construction to param- 
eterize the true initial state with a few parameters which 
can be adjusted to help the rescattering calculation agree 
with the RHIC results. If it is possible to obtain reason- 
able agreement with the data, then one might ask how 
early in the calculation one can go and still maintain a 
physically motivated model. The discussion of this im- 
portant point is deferred until later. 



II. CALCULATIONAL METHOD 

A brief description of the rescattering model calcula- 
tional method is given below. The method used is simi- 
lar to that used in previous calculations for lower CERN 
Super Proton Synchrotron (SPS) energies Rescat- 
tering is simulated with a semi-classical Monte Carlo 
calculation which assumes strong binary collisions be- 
tween hadrons. The Monte Carlo calculation is carried 
out in three stages: 1) initialization and hadronization, 
2) rescattering and freeze out, and 3) calculation of ex- 
perimental observables. Relativistic kinematics is used 
throughout. All calculations are made to simulate RHIC- 
energy Au-|-Au collisions in order to compare with the 
results of the Year-1 RHIC data. 

The hadronization model employs simple parameteri- 
zations to describe the initial momenta and space-time 
of the hadrons similar to that used by Herrmann and 
Bertsch ||l2|. The initial momenta are assumed to follow 
a thermal transverse (perpendicular to the beam direc- 
tion) momentum distribution for all particles, 

{l/mT)dN/dmT = CmT/[exp (mr/T) ± 1] (1) 
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where = \/pt^ + is the transverse mass, pt is 
the transverse momentum, toq is the particle rest mass, C 
is a normalization constant, and T is the initial "temper- 
ature" of the system, and a gaussian rapidity distribution 
for mesons, 



dN/dy = Dexp [-{y - yo) /(2a^')] 



(2) 



where y = Q.^\n[{E p^) / {E — p^)] is the rapidity, E is 
the particle energy, p^ is the longitudinal (along the beam 
direction) momentum, Z? is a normalization constant, 
is the central rapidity value (mid-rapidity), and Uy is the 
rapidity width. Two rapidity distributions for baryons 
have been tried: 1) flat and then falling off near beam ra- 
pidity and 2) peaked at central rapidity and falling off un- 
til beam rapidity. Both baryon distributions give about 
the same results. The initial space-time of the hadrons 
for & = fm (i.e. zero impact parameter or central col- 
lisions) is parameterized as having cylindrical symmetry 
with respect to the beam axis. The transverse particle 
density dependence is assumed to be that of a projected 
uniform sphere of radius equal to the projectile radius, 
R {R — tqA^/^ ^ where ro = 1.12 fm and A is the atomic 
mass number of the projectile). For & > (non-central 
collisions) the transverse particle density is that of over- 
lapping projected spheres whose centers are separated 
by a distance b. The longitudinal particle hadronization 
position {zhad) and time {thad) are determined by the 
relativistic equations p3|, 



Zhad = Thad sinhy; ttad = Thad coshy 



(3) 



where y is the particle rapidity and Thad is the hadroniza- 
tion proper time. Thus, apart from particle multiplici- 
ties, the hadronization model has three free parameters 
to extract from experiment: Uy, T and Thad- The hadrons 
included in the calculation are pions, kaons, nucleons and 
lambdas (tt, K, N, and A), and the p, tj, 77, 77', 0, A, and 
K* resonances. For simplicity, the calculation is isospin 
averaged (e.g. no distinction is made among a 7r+, tt*^, 
and TT^). Resonances are present at hadronization and 
also can be produced as a result of rescattering. Initial 
resonance multiplicity fractions are taken from Herrmann 
and Bertsch p2|, who extracted results from the HELIOS 
experiment The initial resonance fractions used in 
the present calculations are: 77/77 = 0.05, p/ii = 0.1, 
p/io = 3, (j)/{p + w) = 0.12, 77777 = K*/w = 1 and, for 
simphcity, A/A^ = 0. 

The second stage in the calculation is rescattering 
which finishes with the freeze out and decay of all par- 
ticles. Starting from the initial stage (i = fm/c), the 
positions of all particles are allowed to evolve in time in 
small time steps {dt = 0.1 fm/c) according to their ini- 
tial momenta. At each time step each particle is checked 
to see a) if it decays, and b) if it is sufficiently close to 
another particle to scatter with it. Isospin-averaged s- 
wave and p-wave cross sections for meson scattering are 
obtained from Prakash et al. (l5|. The calculation is car- 
ried out to 100 fm/c, although most of the rescattering 
finishes by about 30 fm/c. The rescattering calculation 
is described in more detail elsewhere 
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Calculations are carried out assuming initial parameter 
values and particle multiplicities for each type of parti- 
cle. In the last stage of the calculation, the freeze-out 
and decay momenta and space-times are used to pro- 
duce observables such as pion, kaon, and nucleon mul- 
tiplicities and transverse momentum and rapidity dis- 
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FIG. 1: Pseudorapidity distributions summed over pions, 
kaons, and nucleons from rescattering calculations for 6 = 0, 5, 
and 8 fm. 



tributions. The values of the initial parameters of the 
calculation and multiplicities are constrained to give ob- 
servables which agree with available measured hadronic 
observables. As a cross-check on this, the total kinetic 
energy from the calculation is determined and compared 
with the RHIC center of mass energy of -^i = 130 GeV 
to see that they are in reasonable agreement. Particle 
multiplicities were estimated from the charged hadron 
multiplicity measurements of the RHIC PHOBOS exper- 
iment p^ . Calculations were carried out using isospin- 
summed events containing at freezeout about 5000 pi- 
ons, 500 kaons, and 650 nucleons (A's were decayed). 
The hadronization model parameters used were T = 300 
MeV, ay=2A, and Thad=^ fm/c. It is interesing to note 
that the same value of Thad was required in a previous 
rescattering calculation to successfully describe results 
from SPS Pb-fPb collisions Q. 

Figure 1 shows normalized pseudorapidity distribu- 
tions summed over pions, kaons, and nucleons from 
rescattering calculations for b = 0,5, and 8 fm. The 
widths of the distributions and the flattening near 77 = 
is similar to data from PHOBOS Note that if the 

initial pion rapidity distribution in the calculation was 
slightly flattened at y = from the simple Gaussian in 
Equation 2, the small dip seen in the data at 77 = could 
be reproduced. 

Figure 2 shows ttit distributions for pions, kaons, and 
nucleons from the rescattering calculation for 6 = fm 
near midrapidity (— 1 < y < 1) fitted to exponentials 
of the form exp^—mx/B), where B is the slope param- 
eter. The extracted slope parameters shown in Figure 
2 are close in value to preliminary measurements from 
the STAR experiment for the tt", K~ , and anti-proton 
of 190 ± 10, 300 ± 30, and 565 ± 50 MeV, respectively 
p8[ . Thus, we see that if all hadrons begin at a common 
temperature-parameter value of 300 MeV, the hadronic 
rescattering alone is able to generate enough radial flow 
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FIG. 2: Transverse mass distributions from the rescattering 
model. The lines are exponential fits to the distributions and 
the slope parameters are shown. 

to account for the differences in slope among the pion, 
kaon, and nucleon uit distributions. 

Figure 3 shows transverse momentum distributions for 
pions, kaons, and nucleons from a 6 = 8 fm rescatter- 
ing calculation which extend to high-pr, i.e. 6 GeV/c. 
As observed in PHENIX data , the pion and nucleon 
distributions merge for pT > 2 GeV/c and the kaon dis- 
tribution crosses that for nucleons at around 1 GeV/c. 

The elliptic flow and two-pion HBT observables are 
also calculated from the freeze-out momenta and space- 
time positions of the particles at the end of the rescat- 
tering stage. The elliptic flow variable, V2, is defined as 




V2 = {cos{2(j)))\(j) = arctan(py/pj;) (4) 

where Px and Py are the x and y components of the par- 
ticle momentum, and x is in the impact parameter direc- 
tion and y is in the "out of plane" direction (i.e. a; — z is 
the reaction plane and z is the beam direction). The HBT 
pion source parameters are extracted from the rescatter- 
ing calculation using the same method as was applied for 
previous SPS-energy rescattering calculations The 
Pratt-Bertsch "out-side-long" radius parameterization is 
used ||2l|, ^ yielding the four parameters Rxside, Rrout, 
Rhongi and A, which represent two mutually perpendic- 
ular transverse (to the beam direction) radius parame- 
ters, a radius parameter along the beam direction, and 
a parameter related to the "strength" of the two-pion 
correlations, respectively. 



III. ELLIPTIC FLOW RESULTS 

Figure 4 shows the pt dependence of V2 for pions and 
nucleons extracted from the 6 = 8 fm rescattering cal- 
culation compared with the trends of the STAR mea- 
surements for TT"^ -|- 7r~ and p + p — bar at 11 — 45% 



FIG. 3: Transverse momentum distributions for pions, kaons, 
and nucleons from a 6 = 8 fm rescattering calculation. 

centrality [Q, which roughly corresponds to this impact 
parameter. Figure 5 compares the pt dependence of V2 
for kaons from the 6 = 8 fm rescattering calculation with 
the STAR measurements for s at 11 — 45% central- 
ity As seen, the rescattering calculation values are 
in reasonable agreement with the STAR measurements. 
The flattening out of the pion and nucleon V2 distribu- 
tions for Pt > 2 GeV/c is consistant with that seen in 
STAR and PHENIX results for minimum-bias hadrons 
]l9| , (the kaon V2 calculation does not extend higher 
than 2 GeV/c in Figure 5 due to limited statistics). Thus, 
the same rescattering mechanism that can account for the 
radial flow seen in Figures 2 and 3 is also seen to account 
for the magnitude and pt dependence of the elliptic flow 
for pions, kaons, and nucleons. 

IV. HBT RESULTS 

The pion source parameters extracted from HBT anal- 
yses of rescattering calculations for three different impact 
parameters, 6 = 0, 5, and 8 fm, are compared with STAR 
■K" measurements at three centrality bins j^] in Figure 
6. Note that the PHENIX HBT results |] are in basic 
agreement with the STAR results. The STAR centrality 
bins labeled "3", "2", and "1" in the flgure correspond 
to 12% of central, the next 20%, and the next 40%, re- 
spectively. These bins are roughly approximated by the 
impact parameters used in the rescattering calculations, 
i.e. the average impact parameters of the STAR central- 
ity bins are estimated to be within ±2 fm of the rescatter- 
ing calculation impact parameters used to compare with 
them. In the left panel, the centrality dependence of the 
HBT parameters is plotted for a pt bin of 0.125 — 0.225 
GeV/c. In the right panel, the tot dependence of the 
HBT parameters is plotted for centrality bin 3, for the 
STAR measurments, or 6 = fm, for the rescattering 
calculations. Although there are differences in some of 
the details, the trends of the STAR HBT measurements 
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FIG. 4: Calculations of V2 from the rescattering model for 
6 = 8 fm for pions and nucleons compared with STAR mea- 
surements at 11-45% centrality. The plotted points with error 
bars are the rescattering calculations and the lines show the 
trends of the STAR measurements. Average errors on the 
STAR measurements are < 0.002 for pions and 0.006 for pro- 
tons -l-antiprotons. 



are seen to be described rather well by the rescattering 
calculation. 



V. DISCUSSION 

As shown above, the radial and elliptic flow as well as 
the features of the HBT measurements at RHIC can be 
adequately described by the rescattering model with the 
hadronization model parameters given earlier. The re- 
sults of the calculations are found to be sensitive to the 
value of Thad used, as was studied in detail for SPS rescat- 
tering calculations Q. For calculations with Thad > 1 
fm/c the initial hadron density is smaller, fewer collisions 
occur, and the rescattering-generated flow is reduced, re- 
ducing in magnitude the radial and elliptic flow and most 
of the HBT observables. Only the HBT parameter RLong 
increases for larger r^ad reflecting the increased longitu- 
dinal size of the initial hadron source, as seen in Equa- 
tion 3. One can compensate for this reduced flow in the 
other observables by introducing an ad hoc initial "flow 
velocity parameter", but the increased RLong cannot be 
compensated by this new parameter. In this sense, the 
initial hadron model used in the present calculations with 
Thad ~ 1 fm/c and no initial flow is uniquely determined 

with the help of RLong- 

At this point, one can consider the physical significance 
of the present rescattering model in two different ways. 
The first way is to accept that it is physically valid in the 
time range where hadronic rescattering should be valid, 
e.g. for times later than when the particle density reaches 
about 1 fm~^, and to take the initial state hadroniza- 
tion model as merely a parameterization useful to fit the 
data. Considering the calculation this way, one can at 
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FIG. 5: Calculation of «2 from the rescattering model for 
kaons at b = 8 fm compared with STAR measurements for 
K''s at 11-45% centrality. 



least expect to gain an insight into the phase-space con- 
figuration of the system relatively early in the collision 
(in the calculation 1 fm~^ occurs at a time 4 fm/c after 
hadronization) . The second way to consider the calcula- 
tion is to see if it is possible to also physically motivate 
the initial state hadronization model. An attempt to do 
this is given below. 

In order to consider the present initial state model as 
a physical picture, one must assume: 1) hadronization 
occurs very rapidly after the nuclei have passed through 
each other, i.e. Thad — 1 fm/c, 2) hadrons or at least 
hadron-like objects can exist in the early stage of the 
collision where the maximum value of p approaches 8 
GeV/ fm^, and 3) the initial kinetic energies of hadrons 
can be large enough to be described by T = 300 MeV in 
Equation 1. 

Addressing assumption 2) first, in the calculation the 
maximum number density of hadrons at mid-rapidity at 
t = fm/c is 6.8 fm^^, rapidly dropping to about 1 
fm^^ at i = 4 fm/c. Since most of these hadrons are 
pions, it is useful as a comparison to estimate the effec- 
tive volume of a pion in the context of the tt — tt scat- 
tering cross section, which is about 0.8 frri^ for s-waves 
|]l5| . The "radius" of a pion is found to be 0.25 fm and 
the effective pion volume is 0.065 fm^, the reciprocal of 
which is about 15 fm~^. From this it is seen that at 
the maximum hadron number density in the calculation, 
the particle occupancy of space is estimated to be less 
than 50%, falling rapidly with time. One could specu- 
late that this may be enough spacial separation to allow 
individual hadrons or hadron-like objects to keep their 
identities and not melt into quark matter, resulting in a 
"super- heated" semi-classical gas of hadrons at very early 
times, as assumed in the present calculation. 

Since the calculation takes the point of view of being 
purely hadronic, it is instructive to consider assumptions 
1) and 3) in the context of the Hagedorn thermodynamic 
model of hadronic collisions p5|. According to Hage- 
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FIG. 6: Comparison of HBT source parameters from rescattering with STAR measurements as a function of centrality bin 
(see text) and rriT- The STAR measurements are the solid symbols and the rescattering calculations are the open symbols. 
The errors on the STAR measurements are statistical+systematic. 



dorn, the mass spectrum of hadrons of mass m increases 
proportional to exp(m/ro) in hadronic collisions, where 
To ~ 160 MeV is the limiting temperature of the sys- 
tem. This seems to contradict the value T = 300 MeV 
needed in the present case in Equation 1 to describe the 
data. The Hagedorn model assumes that a) the system 
comes to equilibrium and b) the details of particle pro- 
duction via direct processes and through resonance decay 
average out. Neither of these assumptions is necessarily 
guaranteed at very early times in the collision. The use 
of the thermal functional form, Equation 1 , to set up the 
initial transverse momenta of the hadrons in the present 
calculations is convenient but not required. For example, 
the exponential form exp (— tot/Jc) (where Te is a slope 
parameter) which does not describe thermal equilibrium, 
could have equally well been used. This exponential form 
of the transverse mass distribution was successfully used 
previously in rescattering calculations to describe SPS 
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Assumption 1) can also be motivated by the Color 
Glass Condensate model In the usual version 

of this picture, after the collision takes place the Color 
Glass melts into quarks and gluons in a timescale of about 
0.3 fm/c at RHIC energy, and then the matter expands 
and thermalizes into quark matter by about 1 fm/c. In 
the context of the present rescattering calculations, it is 
tempting to modify the collision senario such that in- 



stead of the Color Glass melting into quarks and gluons 
just after the collision, the sudden impact of the collision 
"shatters" it directly into hadronic fragments on the same 
timescale as in the parton senario due to the hadronic 
strong interactions. 



VI. SUMMARY 

A simple thermal-like initial state model coupled with 
a hadronic rescattering calculation is able to adequately 
describe the large elliptic flow and small HBT radii re- 
cently measured at RHIC. A feature of this picture is a 
very early hadronization time of about 1 fm/c after the 
initial collision of the nuclei. 
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